EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 



Ch- LHCb-PAPER-201 1-027 

CERN-PH-EP-2012-039 
February 29, 2012 



(N 

O 

(N 

^ Opposite-side flavour tagging of B 
S mesons at the LHCb experiment 

The LHCb collaboration [3 

(N 

^ Abstract 

The cahbration and performance of the opposite-side flavour tagging algorithms 
0^ used for the measurements of time-dependent asymmetries at the LHCb experiment 

are described. The algorithms have been developed using simulated events and op- 
g timized and calibrated with B+ J/iIjK+, J/ipK*° and 5° D*-^i+v^ 

decay modes with 0.37 fb~^ of data collected in pp collisions at -^/s = 7TeV dur- 
I ing the 2011 physics run. The opposite-side tagging power is determined in the 

j> — > J/ipK^ channel to be (2.10±0.08±0.24)%, where the first uncertainty is 
statistical and the second is systematic. 



Submitted to Eur. Phys. J. C 



Authors arc listed on the following pages. 



ii 



The LHCb collaboration 



R. Aaij38, C. Abellan Beteta^^-", B. Adeva3^ M. Adinolfi^^^ C. Adrover^ A. AfFolder^'', 
Z. Ajaltouni^, J. Albrecht^^, F. Alessio^^, M. Alexander"'^, G. Alkhazov^'', 

P. Alvarez Cartelle^*^, A.A. Alves Jr^^, S. Amato^, Y. Amhis^^, J. Anderson^'', R.B. Appleby^\ 
O. Aquines Gutierrez^'^, F. Archilli-^^'^^, L. Arrabito^^, A. Artamonov M. Artuso^^'^^, 

E. Aslanides^, G. Auriemma^^'™, S. Bachmann^^, J.J. Back"'^, D.S. Bailey^^, V. Balagura^*'^^ , 
W. Baldini^^ R.J. Barlow^^, C. Barschcl'■^■^ S. Barsuk^ W. Barter''^, A. Bates^^ C. Baucr^o, 
Th. Bauer^*^, A. Bay^*^, I. Bediaga^ S. Belogurov^*^, K. Belous^^, I. Belyaev^^, E. Ben-Haim®, 
M. Benayoun^, G. Bencivenni^^, S. Benson^^, J. Benton^^, R. Bernet^^, M.-O. Bettler^^, 

M. van Beuzekom^s, A. Bien^^, S. Bifani^^ T. Bird^^, A. Bizzeti^^''', P.M. Bj0rnstad5^ 

T. Blakc^s, F. Blanc^^ C. Blanks^", J. Blouw", S. Blusk^^ A. Bobrov^i, V. Bocci^^, 

A. Bondar^^ N. Bondar^'', W. Bonivento^^, S. Borghi'^^'^^ A. Borgia^^, T.J.V. Bowcock"'^, 

C. Bozzi^^, T. Brambach^, J. van den Brand^'^, J. Bressieux^^, D. Brett^S M. Britsch^°, 

T. Britton^'^, N.H. Brook^^, H. Brown^^, K. dc Bruyn'^^, A. Biichler-Germann^''', I. Burducea^^, 

A. Bursche-'^'^, J. Buytaert^^, S. Cadeddul^ O. Callot^ M. Calvi^^'^ M. Calvo Gomez^^.", 

A. Camboni33, R Campana^^-^s , A. Carbone^"^, G. Carboni^^''^, R. Cardinalel9'^'^^ 

A. Cardiml^ L. Carson^^, K. Carvalho Akiba^, G. Casse^^^ ^ Cattaneo3^ Ch. Cauet^ 

M. Charlcs^^ Ph. Charpcntier3^ N. Chiapolinr'^^ K. Ciba^^, X. Cid VidaP, G. Ciczarek^o, 

P.E.L. Clarkc^'^'^s, M. Clemencic'^-'\ H.V. Cliff'^'^, J. Closier^s, C. Coca^^, V. Cocoas, J. Cogan^, 

P. Collins^^, A. Comerma-Montells^^, F. Constantin^*^, A. Contu^^, A. Cook*^^, M. Coombes^^, 

G. Corti^s, B. Couturier^s^ q ^ Cowan^^, R. Currie^^, C. D'Ambrosio^s, P. David^, 
P.N.Y. David^s, I. De Bonis^ S. De Capua^i''^, M. Dc Cian'^^ F. Dc Lorcnzi^^ 
J.M. Dc Miranda^, L. Dc Paula^, P. Dc Simonc^^, D. Decamp-^, M. Dcckcnhoff, 

H. Degaudenzi^^'^^, L. Del Buono^, C. Deplano^^, D. Derkacli^'*''^^ , O. Deschamps^, 

F. Dettori^*^, J. Dickens^^, H. Dijkstra^^, P. D iniz Batista^, F. Domingo Bonal^^'", 

S. Donleavy'^^, F. Dordei^\ A. Dosil Suarez^'', D. Dossctt'^^ A. Dovbnya^°, F. Dupcrtuis^^, 
R. Dzhelyadin^^, A. Dziurda^^, S. Easo"'^, U. Egede^°, V. Egorychev^*, S. Eidelman^^ 

D. van Eijk^s, F. Eisele", S. Eisenhardt^^, R. Ekelhof^, L. Eklund^^^, Ch. Elsasser^^, 
D. Elsby^2^ D Esperante Pereira^^, A. Falabella^^'^'^^, E. Fanchini^^'-' , C. Farber^\ 

G. Fardell^^, C. Farinelli^^, S. Farry^^^ y pave^*^, V. Fernandez Albor^^, M. Ferro-Luzzi^^, 
S. Fihppov^o, C. Fitzpatrick^^, M. Fontana^°, F. FontaneUi^^-*, R. Forty3■^ O. Francisco^, 
M. Frank^s, C. Frei^^, M. Frosini^'^'-^, S. Furcas^o, A. Gallas Torreira^^, D. Gam^^'^, 

M. Gandelman^, P. Gandini^^^ Y. Gao=^, J-C. Garnier^^ J. Garofoli^^, J. Garra Tico''^, 

L. Garrido^^ D. Gascon^^, C. Gaspares, R. Gauld^^ N. Gauvin3^ M. Gersabeck^^, 

T. Gershon^5'35, Ph. Ghez^, V. Gibson^^, V.V. GHgorov^s, C. G6bel54, jy Golubkov^^, 

A. Golutvin^°'^*'^^, A. Gomes^, H. Gordon^^, M. Grabalosa Gandara^^, R. Graciani Diaz^^, 

L.A. Granado Cardoso^^, E. Grauges^^, G. Graziani^'^, A. Grecu^^, E. Greening^^, 

S. Gregson^^, B. Gui^^, E. Gushchin^o, Yu. Guz^^, T. Gys^^, C. Hadjivasiliou^^ G. Haefeli3^ 

C. Haen^^, S.C. Haines^^, T. Hampson^^, S. Hansmann-Mcnzemer^^, R. Harji^*^, N. Harnew^^, 
J. Harrison^^, P.F. Harrison"^^, T. Hartmann^^, J. He^, V. Heijne^^, K. Hennessy^^, 

P. Henrard^, J.A. Hernando Morata^^, E. van Herwijnen^^, E. Hicks^^, K. Holubyev^^, 
P. Hopchev^, W. Hulsbergen3^ P. Hunt^^^ T. Huse^^, R.S. Huston^^ D. Hutchcroft^^ 

D. Hynds^^, V. lakovenko^^, P. Ilten^^, J. Imong^^, R. Jacobsson^^, A. Jaeger^^, 

M. Jahjah Hussein^, E. Jans^*, F. Jansen^^, P. Jaton^'^, B. Jean-Marie'^, F. Jing^, M. John^^, 
D. Johnson52, C.R. Jones^^, B. Jost^^, M. Kaballo^ S. Kandybei^o, M. Karacson^s, 



iii 



T.M. Karbach^, J. Keaveney^^ I.R. Kenyon^^ U. Kerzel^^, T. Ketel^^, A. Keune^^, 

B. Khanji^, Y.M. Kim^^, M. Knecht^^, R.F. Koopman^'', P. Koppenburg^^, M. Korolev^^, 

A. Kozlinskiy^^, L. Kravchuk^'', K. Kreplin^^, M. Kreps^^, G. Krocker^^, P. Krokovny^^, 

F. Kruse^ K. Kruzelecki^^ , M. Kucharczyk^O'^s.ssj^ t_ Kvaratskheliya^^'^s , V.N. La Thi^^, 
D. Lacarrere^^, G. Lafferty^\ A. Lai^^, D. Lambert"'^, R.W. Lambert"^^, E. Lanciotti^^, 

G. Lanfranchi^^, C. Langenbruch^^, T. Latham^^, C. Lazzeroni^^, R. Le Gac^, 

J. van Lccrdam'^^, J. -P. Lccs'^, R. Lefevre^, A. Leflat'^^'^^, J. Lefrancois''', O. Leroy^, 

T. Lesiak23, L. Li^, L. Li Gioi^ M. Licng^, M. Lilcs'^^, R. Lindner^s, C. Linn", B. Liu^ 

G. Liu^^, J. von Loeben^°, J.H. Lopes^, E. Lopez Asamar"^^, N. Lopez-March^®, H. Lu^, 

J. Luisier^^, A. Mac Raighne^^, F. Machefert''', LV. Machikhiliyan'^'^s, F. Maciuc^°, 

O. Maev27'35, J. Magnin\ S. Maldc'^^ R.M.D. Mamunur3^ G. Mancai^''^, G. Mancinelli®, 

N. Mangiafave^^, U. Marconi^'i, R. Marki^*^, J. Marksi\ G. Martellotti22, A. Martens^, 

L. Martin^^, A. Martm Sanchez'^, D. Martinez Santos^^, A. Massafferri^ Z. Mathe^^, 

C. Matteuzzi^o, M. Matveev^''', E. Maurice®, B. Maynard^^, A. Mazurov^®'^^'^^, G. McGregor^^ 
R. McNulty^^, M. Meissner", M. Merk^^, J. Merkel^, R. Messi^^'^, S. Miglioranzi^^, 

D. A. Milanes^^, M.-N. Minard"^, J. Molina Rodriguez^^, S. Montcil^, D. Moran^^, 

P. Morawski23, R. Mountain^^, L Mous^^, F. Muheim'^'^, K. Miiller^^, R. Muresan^®, 

B. Muryn24, B. Muster^®, M. Musy^^, J. Mylroie-Smith^^, P. Naik4^ T. Nakada^®, 

R. Nandakumar^®, L Nasteva^ M. Nedos^, M. Needham''^, N. Neufeld^^, A.D. Nguyen^®, 

C. Nguyen- Mau3®'°, M. Nicof , V. Nicss^, N. Nikitin^'^, A. Nomerotski^2,35^ Novoselov^^, 
A. Oblakowska-Mucha^^, V. Obraztsov^^ S. O ggero^^, S. Ogilvy^®, O. Okhrimenko^^ , 

R. 01demani5,d,35^ Orlandea^®, J.M. Otalora Goicochea^, P. Owen^^, K. Pal5^ 

J. Palacios^''', A. Palano^^''', M. Palutan^^, J. Panman^^, A. Papanestis^®, M. Pappagallo^^, 

C. Parkes5\ C.J. Parkinson^o, G. Passaleva^'^, G.D. Patel^^ M. PateP, S.K. Paterson^^, 
G.N. Patrick^®, C. Patrignani^^'*, C. Pavel-Nicorescu^®, A. Pazos Alvarez^^, A. Pellegrino^^, 
G. Penso^^'', M. Pepe Altarelli^^, S. Perazzini^^'^ D.L. Perego^^'-?, E. Perez Trigo^^, 

A. Perez-Calero Yzquierdo^^, P. Perret^, M. Perrin-Terrin®, G. Pessina^°, A. Petrella^®'^^, 

A. Petrolinii9'«, A. Phan^^^ Picatoste Olloqui^^, B. Pie Valls^^, B. Pietrzyk"^, T. Pilaf^^, 

D. Pinci22, R. Plackett'^^ S. Playfer'^^ M. Plo Casasus^^, G. Polok^^, A. Poluektov'^^-^i, 

E. Polycarpo^, D. Popov^°, B. Popovici^®, C. Potterat^^, A. Powell^^, J. Prisciandaro^®, 

V. Pugatch^^, A. Puig Navarro^^, W. Qian^^, J.H. Rademacker^^, B. Rakotomiaramanana^®, 

M.S. Rangcl2, L Raniuk^^ G. Ravcn^^, S. Redford52, j^i.m. Reid^^^^ A.C. dos Reis^ 

S. Ricciardi-*®, A. Richards^^^, K. Rinnert'^^ D.A. Roa Romero^ P. Robbe^ E. Rodrigues^^^'^S 

F. Rodrigues^, P. Rodriguez Perez^^, G.J. Rogers''^, S. Roiser^^, V. Romanovsky^^ , 

M. Rosello^^.n^ J Rouvinet^®, T. Ruf^, H. Ruiz^^, G. Sabatino^^'*^, J.J. Saborido Silva3^ 

N. Sagidova^^, P. Sail^^, B. Saitta^^'^^, C. Salzmann^''', M. Sannino^^'S R. Santacesaria^^^ 

C. Santamarina Rios^^, R. Santinelli2^ E. Santovetti^^'^, M. Sapunov®, A. Sarti^^-', 

C. Satriano22."^, A. Satta^i, M. Savriel®'^ D. Savrina^^, P. Schaack^^, M. Schiller^s, 

S. Schleich^, M. Schlupp^, M. Schmellingio, B. Schmidt3^ O. Schneider^®, A. Schopper^s, 

M.-H. Schune^, R. Schwemmer^^ , B. Sciascia"*^^, A. Sciubba^S'', M. Seco^ , A. Semennikov^^, 

K. Senderowska^"^, L Sepp^®, N. Serra^^, J. Serrano®, P. Seyfert", M. Shapkin^^, 

L Shapoval^O'^s^ P. Shatalov^^, Y. Shcheglov^^, T. Shears^^ L. Shekhtman^S O. Shevchenko^^, 

V. Shevchenko^s, A. Shires^®, R. Silva Coutinho^^ T. Skwarnicki^a, N.A. Smith^^ 

E. Smith52,46^ Sobczak^ F.J.P. Solcr4^ A. Solomin4^ F. Soomro^^'-'^^ B. Souza De Paula^, 

B. Spaan^, A. Sparkes^^ P. Spradlin4^ F. Stagni^s, S. Stahl", O. Steinkamp^^, S. Stoica^®, 
S. Stone53,35^ B Storaci^^, M. Straticiuc^®, U. Straumann^^, V.K. Subbiah^^, S. Swientek^, 



iv 



M. Szczekowski^s, p. Szczypka^^, T. Szumlak^^, S. T'Jampens"^, E. Teodorescu^^, F. Teubert^^, 

C. Thomas^^, E. Thomas^^ J. van Tilburg^\ V. Tisserand^, M. Tobin^'^, S. Topp-Joergensen^^, 
N. Torr^^, E. Tourncfier'^'^*', S. Tourneur^^, M.T. Tran^®, A. Tsarcgorodtsev^, N. Tuning^^, 

M. Ubeda Garcia^^, A. Ukleja^^, p. Urquijo^^ U. Uwer", V. Vagnomi"^, G. Valenti^^ 
R. Vazquez Gomez^^, P. Vazquez Regueiro^*^, S. Vecchi^^, J.J. Velthuis*^^, M. Veltri^^'^', 
B. Viaud''', I. Videau''', D. Vieira^, X. Vilasis-Cardona^^'", J. Visniakov^^, A. VoUhardt^''', 

D. Volyanskyy^°, D. Voong^^, A. Vorobyev^''', H. Voss^'^, S. Wandernoth^^, J. Wang^''^, 

D.R. Ward4^ N.K. Watson^^ A.D. Webber^i, D. Websdale^o, M. Whitehead''^ , D. Wicdncrii, 
L. Wiggers^s, G. Wilkinson^^, M.P. Williams'^^'^^ M. Williams^o, F.F. Wilson-^^ J. Wishahi^ 
M. Witek23, W. Witzeling^s, S.A. Wotton^^, K. Wyllie^s, Y. Xie^^, F. Xing^^, Z. Xing^^, 
Z. Yang3, R. Young'^'^, O. Yushchenko^^ , M. Zangoli^^, M. Zavcrtyacv^°'", F. Zhang^, 
L. Zhang^^, W.C. Zhang^^, Y. Zhang^, A. Zhelezov^^ L. Zhong^, A. Zvyagin^^. 

^ Centra Brasileiro de Pesquisas Fisicas ( CBPF), Rio de Janeiro, Brazil 
Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil 
^ Center for High Energy Physics, Tsinghua University, Beijing, China 
^LAPP, Universite de Savoie, CNRS/IN2P3, Annecy-Le-Vieux, France 

^Clermont Universite, Universite Blaise Pascal, CNRS/IN2P3, LPC, Clermont-Ferrand, France 
^CPPM, Aix-Marseille Universite, CNRS/IN2P3, Marseille, France 

"^LAL, Universite Paris-Sud, CNRS/IN2P3, Orsay, France 

^LPNHE, Universite Pierre et Marie Curie, Universite Paris Diderot, CNRS/IN2P3, Paris, France 

^Fakultdt Physik, Technische Universitdt Dortmund, Dortmund, Germany 
Max- Planck- Institut fiir Kernphysik (MPIK), Heidelberg, Germany 
Physikalisches Institut, Ruprecht-Karls- Universitdt Heidelberg, Heidelberg, Germany 
School of Physics, University College Dublin, Dublin, Ireland 

^^Sezione INFN di Bari, Bari, Italy 

^'^Sezione INFN di Bologna, Bologna, Italy 

^^Sezione INFN di Cagliari, Cagliari, Italy 

^^Sezione INFN di Ferrara, Ferrara, Italy 

^'^Sezione INFN di Firenze, Firenze, Italy 

^^Laboratori Nazionali dell'INFN di Frascati, Frascati, Italy 

^^Sezione INFN di Genova, Geneva, Italy 

"^^Sezione INFN di Milano Bicocca, Milano, Italy 

"^^Sezione INFN di Roma Tor Vergata, Roma, Italy 

"^^Sezione INFN di Roma La Sapienza, Rom,a, Italy 

'^^Henryk Niewodniczanski Institute of Nuclear Physics Polish Academy of Sciences, Krakow, Poland 
'^'^AGH University of Science and Technology, Krakow, Poland 

'^^ Saltan Institute far Nuclear Studies, Warsaw, Poland 

^^Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest- Magurele, Romania 

'^'^ Petersburg Nuclear Physics Institute (PNPI), Gatchina, Russia 

'^^ Institute of Theoretical and Experimental Physics (ITEP), Moscow, Russia 

'^^ Institute of Nuclear Physics, Moscow State University (SINP MSU), Moscow, Russia 

Institute for Nuclear Research of the Russian Academy of Sciences (INR RAN), Moscow, Russia 
^^Budker Institute of Nuclear Physics (SB RAS) and Novosibirsk State University, Novosibirsk, Russia 
^"^ Institute for High Energy Physics (IHEP), Protvino, Russia 

Universitat de Barcelona, Barcelona, Spain 

Universidad de Santiago de Compostela, Santiago de Compostela, Spain 
European Organization for Nuclear Research (CERN), Geneva, Switzerland 
^^Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne, Switzerland 

^"^ Physik- Institut, Universitdt Ziirich, Zurich, Switzerland 

^^Nikhef National Institute for Subatomic Physics, Amsterdam, The Netherlands 



V 



^^Nikhef National Institute for Subatomic Physics and Vrije Universiteit, Amsterdam, The Netherlands 
^^NSC Kharkiv Institute of Physics and Technology (NSC KIPT), Kharkiv, Ukraine 

Institute for Nuclear Research of the National Academy of Sciences (KINR), Kyiv, Ukraine 

University of Birmingham, Birmingham, United Kingdom 
^^H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom 
^'^ Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom 

Department of Physics, University of Warwick, Coventry, United Kingdom 
^^STFC Rutherford Appleton Laboratory, Didcot, United Kingdom 

'^'^ School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom 
'^^ School of Physics and Astronomy, University of Glasgow, Clasgow, United Kingdom 

Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom 
^'^ Imperial College London, London, United Kingdom 

School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom 

Department of Physics, University of Oxford, Oxford, United Kingdom 

Syracuse University, Syracuse, NY, United States 
^^Pontiftcia Universidade Catolica do Rio de Janeiro (PUC-Rio), Rio de Janeiro, Brazil, associated to ^ 

CC-IN2P3, CNRS/IN2P3, Lyon- Villeurbanne, France, associated member 

Physikalisches Institut, Universitdt Rostock, Rostock, Germany, associated to 

"■P.N. Lebedev Physical Institute, Russian Academy of Science (LPI RAS), Moscow, Russia 

^ Universitd di Bari, Bari, Italy 

Universita di Bologna, Bologna, Italy 

Universitd di Cagliari, Cagliari, Italy 

Universita di Ferrara, Ferrara, Italy 
^ Universita di Firenze, Firenze, Italy 
^Universita di Urbino, Urbino, Italy 
^ Universita di Modena e Reggie Emilia, Modena, Italy 
* Universita di Genova, Genova, Italy 
^ Universita di Milano Bicocca, Milano, Italy 
^ Universita di Roma Tor Vergata, Roma, Italy 
' Universita di Roma La Sapienza, Roma, Italy 

Universita della Basilicata, Potenza, Italy 
"LIFAELS, La Salle, Universitat Ramon Hull, Barcelona, Spain 
"Hanoi University of Science, Hanoi, Viet Nam 



vi 



1 Introduction 



The identification of tlie flavour of reconstructed and mesons at production is 
necessary for tlie measurements of oscillations and time-dependent CP asymmetries. This 
procedure is known as flavour tagging and is performed at LHCb by means of several 
algorithms. 

Opposite-side (OS) tagging algorithms rely on the pair production of b and b quarks 
and infer the flavour of a given B meson (signal B) from the identification of the flavour 
of the other b hadroi][^ (tagging B). The algorithms use the charge of the lepton (/i, e) 
from semileptonic b decays, the charge of the kaon from the 6 — >■ c — )■ s decay chain or the 
charge of the inclusive secondary vertex reconstructed from 6-hadron decay products. All 
these methods have an intrinsic dilution on the tagging decision, for example due to the 
possibility of flavour oscillations of the tagging B. This paper describes the optimization 
and calibration of the OS tagging algorithms which are performed with the data used 
for the first measurements performed by LHCb on B^ mixing and time-dependent CP 
violation (l}]3]. 

Additional tagging power can be derived from same-side tagging algorithms which 
determine the flavour of the signal B by exploiting its correlation with particles produced 
in the hadronization process. The use of these algorithms at LHCb will be described in 
a forthcoming publication. The use of flavour tagging in previous experiments at hadron 
colliders is described in Refs. (4|[5]. 

The sensitivity of a measured CP asymmetry is directly related to the effective tagging 
efficiency Ees, or tagging power. The tagging power represents the effective statistical 
reduction of the sample size, and is defined as 

£eS = ^tag'P^ = £tag(l " 2c<;)^ (1) 

where Stag is the tagging efficiency, u is the mistag fraction and V is the dilution. The 
tagging efficiency and the mistag fraction are defined as 

R + W W 

where R, W, U are the number of correctly tagged, incorrectly tagged and untagged 
events, respectively. 

The mistag fraction can be measured in data using flavour-specific decay channels, i.e. 
those decays where the final state particles uniquely define the quark/antiquark content 
of the signal B. In this paper, the decay channels B^ — )■ J/ipK^, B^ — )■ J/ipK*^ and 
5"— 7- D*^/i+z/^ are used. For charged mesons, the mistag fraction is obtained by directly 
comparing the tagging decision with the flavour of the signal S, while for neutral mesons 
it is obtained by fitting the B^ flavour oscillation as a function of the decay time. 

The probability of a given tag decision to be correct is estimated from the kinematic 
properties of the tagging particle and the event itself by means of a neural network 



^Unless explicitly stated, charge conjugate modes are always included throughout this paper. 
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trained on Monte Carlo (MC) simulated events to identify the correct flavour of the 
signal B. When more than one tagging algorithm gives a response for an event, the 
probabilities provided by each algorithm are combined into a single probability and the 
decisions are combined into a single decision. The combined probability can be exploited 
on an event-by-event basis to assign larger weights to events with low mistag probability 
and thus to increase the overall significance of an asymmetry measurement. In order 
to get the best combination and a reliable estimate of the event weight, the calculated 
probabilities are calibrated on data. The default calibration parameters are extracted 
from the — t- J/ipK^ channel. The other two flavour-specific channels are used to 
perform independent checks of the calibration procedure. 

2 The LHCb detector and the data sample 

The LHCb detector j6] is a single-arm forward spectrometer which measures CP violation 
and rare decays of hadrons containing b and c quarks. A vertex detector (VELO) deter- 
mines with high precision the positions of the primary and secondary vertices as well as 
the impact parameter (IP) of the reconstructed tracks with respect to the primary ver- 
tex. The tracking system also includes a silicon strip detector located in front of a dipole 
magnet with integrated fleld about 4 Tm, and a combination of silicon strip detectors and 
straw drift chambers placed behind the magnet. Charged hadron identification is achieved 
through two ring-imaging Cherenkov (RICH) detectors. The calorimeter system consists 
of a preshower detector, a scintillator pad detector, an electromagnetic calorimeter and 
a hadronic calorimeter. It identifies high transverse energy hadron, electron and photon 
candidates and provides information for the trigger. Five muon stations composed of 
multi-wire proportional chambers and triple-GEMs (gas electron multipliers) provide fast 
information for the trigger and muon identification capability. 

The LHCb trigger consists of two levels. The first, hardware-based, level selects lep- 
tons and hadrons with high transverse momentum, using the calorimeters and the muon 
detectors. The hardware trigger is followed by a software High Level Trigger (HLT), 
subdivided into two stages that use the information from all parts of the detector. The 
first stage performs a partial reconstruction of the event, reducing the rate further and 
allowing the next stage to fully reconstruct and to select the events for storage up to a 
rate of 3 kHz [T] . 

The majority of the events considered in this paper were triggered by a single hadron 
or muon track with large momentum, transverse momentum and IP. In the HLT, the 
channels with a J/tp meson in the final state were selected by a dedicated di-muon decision 
that does not apply any requirement on the IP of the muons. 

The data used in this paper were taken between March and June 2011 and correspond 
to an integrated luminosity of 0.37 fb~ . The polarity of the LHCb magnet was reversed 
several times during the data taking period in order to minimize systematic biases due to 
possible detector asymmetries. 
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3 Flavour tagging algorithms 



Opposite-side tagging uses the identification of electrons, muons or kaons that are at- 
tributed to the other b hadron in the event. It also uses the charge of tracks consistent 
with coming from a secondary vertex not associated with either the primary or the signal B 
vertex. These taggers are called electron, muon, kaon and vertex charge taggers, respec- 
tively. The tagging algorithms were developed and studied using simulated events [s]. 
Subsequently, the criteria to select the tagging particles and to reconstruct the vertex 
charge are re-tuned, using the J/ipK^ and the 5°—)- D*"//"*"//^ control channels. 

An iterative procedure is used to find the selection criteria which maximize the tagging 
power eeff- 

Only charged particles reconstructed with a good quality of the track fit are used. In 
order to reject poorly reconstructed tracks, the track is required to have a polar angle 
with respect to the beamline larger than 12 mrad and a momentum larger than 2 GeV/ c. 
Moreover, in order to avoid possible duplications of the signal tracks, the selected particles 
are required to be outside a cone of 5 mrad formed around any daughter of the signal B. 
To reject tracks coming from other primary interactions in the same bunch crossing, the 
impact parameter significance with respect to these pile-up (PU) vertices, IPpu/o"iPpu > 3, 
is required. 

3.1 Single-particle taggers 

The tagging particles are selected exploiting the properties of the 6-hadron decay. A large 
impact parameter significance with respect to the primary vertex (IP/cxip) and a large 
transverse momentum are required. Furthermore, particle identification cuts are used 
to define each tagger based on the information from the RICH, calorimeter and muon 
systems. For this purpose, the differences between the logarithm of the likelihood for the 
muon, electron, kaon or proton and the pion hypotheses (referred as DLL^_7r, DLLe_7r, 
DLLa'-tt and DLLp_7r) are used. The detailed list of selection criteria is reported in Ta- 
ble [Tj Additional criteria are used to identify the leptons. Muons are required not to share 
hits in the muon chambers with other tracks, in order to avoid mis-identification of tracks 
which are close to the real muon. Electrons are required to be below a certain threshold 
in the ionization charge deposited in the silicon layers of the VELO, in order to reduce the 
number of candidates coming from photon conversions close to the interaction point. An 
additional cut on the ratio of the particle energy E as measured in the electromagnetic 
calorimeter and the momentum p of the candidate electron measured with the tracking 
system, E /p > 0.6, is applied. 

In the case of multiple candidates from the same tagging algorithm, the single-particle 
tagger with the highest is chosen and its charge is used to define the fiavour of the 
signal B. 
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Table 1: Selection criteria for the OS muon, electron and kaon taggers. 



Tagger 


min pt 

[ GeV/c ] 


min p 

[GeV/c] 


min (IP/(Tip) 


Particle identification 
cuts 


min (IPpu/(TiPpu) 


fi 


1.2 


2.0 




DLL^_^ > 2.5 


3.0 


e 


1.0 


2.0 


2.0 


DLLe-^ > 4.0 


3.0 


K 


0.8 


5.9 


4.0 


T)\A.K-^ > 6.5 
DLLk-p > -3.5 


4.7 



3.2 Vertex charge tagger 

The vertex charge tagger is based on the inclusive reconstruction of a secondary vertex 
corresponding to the decay of the tagging B. The vertex reconstruction consists of build- 
ing a composite candidate from two tracks with a transverse momentum > 0.15 GeV/c 
and IP/cTip > 2.5. The pion mass is attributed to the tracks. Moreover, good quality of 
the vertex reconstruction is required and track pairs with an invariant mass compatible 
with a Kg meson are excluded. For each reconstructed candidate the probability that it 
originates from a 6-hadron decay is estimated from the quality of the vertex fit as well 
as from the geometric and kinematic properties. Among the possible candidates the one 
with the highest probability is used. Tracks that are compatible with coming from the 
two track vertex but do not originate from the primary vertex are added to form the final 
candidate. Additional requirements are applied to the tracks asspociated to the recon- 
structed secondary vertex: total momentum > 10 GeV/c, total pt > 1-5 GeV/c, total 
invariant mass > 0.5 GeV/c^ and the sum of IP/crip of all tracks > 10. 

Finally, the charge of the tagging B is calculated as the sum of the charges Qi of all the 
tracks associated to the vertex, weighted with their transverse momentum to the power k 

Qv. = (3) 

where the value n = 0.4 optimizes the tagging power. Events with IQvtxl < 0.275 are 
rejected as untagged. 

3.3 Mistag probabilities and combination of taggers 

For each tagger i, the probability rji of the tag decision to be wrong is estimated by 
using properties of the tagger and of the event itself. This mistag probability is evaluated 
by means of a neural network trained on simulated B~^ — )■ J/ip events to identify the 
correct flavour of the signal B and subsequently calibrated on data as explained in Sect. |5] 
The inputs to each of the neural networks are the signal B transverse momentum, the 
number of pile-up vertices, the number of tracks preselected as tagging candidates and 
various geometrical and kinematic properties of the tagging particle (p, px and IP/o"ip of 
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the particle), or of the tracks associated to the secondary vertex (the average values of 
PT) of IP, the reconstructed invariant mass and the absolute value of the vertex charge). 

If there is more than one tagger available per event, the decisions provided by all 
available taggers are combined into a final decision on the initial flavour of the signal B. 
The combined probability P{b) that the meson contains a 6-quark is calculated as 

where 

Pib) = n - Ml - V^?) , Pib) = n + - ■ (5) 

Here, di is the decision taken by the i-th tagger based on the charge of the particle with 
the convention di = 1(— 1) for the signal B containing a b{b) quark and r/i the correspond- 
ing predicted mistag probability. The combined tagging decision and the corresponding 
mistag probability are d = —1 and r] = 1 — P{h) if P{h) > P(b), otherwise d = +1 and 
ri= 1-P{b). 

The contribution of taggers with a poor tagging power is limited by requiring the 
mistag probabilities of the kaon and the vertex charge to be less than 0.46. 

Due to the correlation among taggers, which is neglected in Eq. [5} the combined 
probability is slightly overestimated. The largest correlation occurs between the vertex 
charge tagger and the other OS taggers, since the secondary vertex may include one 
of these particles. To correct for this overestimation, the combined OS probability is 
calibrated on data, as described in Sect. |5| 



4 Control channels 

The flavour-specific B decay modes 5+ — J/ipK^, B^ — )■ J/tpK*^ and B^ — )■ are 
used for the tagging analysis. All three channels are useful to optimize the performance of 
the OS tagging algorithm and to calibrate the mistag probability. The first two channels 
are chosen as representative control channels for the decays -B°— >■ J/ip4' cind -8°—)- J/ip fo, 
which are used for the measurement of the B^ mixing phase (ps p||3] , and the last channel 
allows detailed studies given the high event yield of the semileptonic decay mode. All 
B decay modes with a J/tp meson in the final state share the same trigger selection 
and common offline selection criteria, which ensures a similar performance of the tagging 
algorithms. Two trigger selections are considered, with or without requirements on the IP 
of the tracks. They are labelled "lifetime biased" and "lifetime unbiased" respectively. 



4.1 Analysis of the S+ J/il^K'^ channel 

The B^ — )■ J/ipK^ candidates are selected by combining J/ip ^ and candi- 

dates. The J/ip mesons are selected by combining two muons with transverse momenta 
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Pt > 0.5 GeV/c that form a common vertex of good quality and have an invariant mass 
in the range 3030 — 3150 MeV/ c^. The candidates are required to have transverse mo- 
menta Pt > 1 GeV/ c and momenta p > 10 GeV/ c and to form a common vertex of good 
quahty with the J/ ip candidate with a resulting invariant mass in a window ±90 MeV/ 
around the mass. Additional requirements on the particle identification of muons 
and kaons are applied to suppress the background contamination. To enhance the sample 
of signal events and reduce the dominant background contamination from prompt J/ip 
mesons combined with random kaons, only the events with a reconstructed decay time of 
the candidate t > 0.3 ps are selected. The decay time t and the invariant mass m of 
the B~^ meson are extracted from a vertex fit that includes a constraint on the associated 
primary vertex, and a constraint on the J/ijj mass for the evaluation of the J/ipK invariant 
mass. In case of multiple B candidates per event, only the one with the smallest vertex 
fit is considered. 

The signal events are statistically disentangled from the background, which is domi- 
nated by partially reconstructed 6-hadron decays to J/ipK^X (where X represents any 
other particle in the decay), by means of an unbinned maximum likelihood fit to the 
reconstructed B^ mass and decay time. In total ~ 85 000 signal events are selected with 
a background to signal ratio BjS ^ 0.035, calculated in a window of ±40 MeV/c^ cen- 
tred around the B^ mass. The mass fit model is based on a double Gaussian distribution 
peaking at the B^ mass for the signal and an exponential distribution for the background. 
The time distributions of both the signal and the background are assumed to be exponen- 
tial, with separate decay constants. The fraction of right, wrong or untagged events in 
the sample is determined according to a probability density function (PDF), P(r), that 
depends on the tagging response r, defined by 

{£^tag(l — ^) r= "right tag decision" 
Stag oj r= "wrong tag decision" (6) 

1 ~ £^tag r="no tag decision". 

The parameters oj and £tag (defined in Eq. |2]) are different for signal and background. 
Fig. [l] shows the mass distribution of the selected and tagged events, together with the 
superimposed fit. 

4.2 Analysis of the )- channel 

The B^ — !■ D*^fi^i'^ channel is selected by requiring that a muon and the decay 
D*~ — i- D^{^ K^'K^)'K~ originate from a common vertex, displaced with respect to 
the pp interaction point. The muon and transverse momenta are required to be larger 
than 0.8 GeV/c and 1.8 GeV/c respectively. The selection criteria exploit the long B^ 
and lifetimes by applying cuts on the impact parameters of the daughter tracks, on 
the pointing of the reconstructed B^ momentum to the primary vertex, on the difference 
between the z coordinate of the B^ and vertices, and on the flight distance. Addi- 
tional cuts are applied on the muon and kaon particle identiflcation and on the quality of 
the flts of all tracks and vertices. In case of multiple B candidates per event the one with 
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Figure 1: Mass distribution of OS tagged 5+ — )- J/ipK^ events. Black points are data, the 
solid blue line, red dotted line and green area are the overall fit, the signal and the background 
components, respectively. 

the smallest impact parameter significance with respect to the primary vertex is consid- 
ered. Only events triggered in the HLT by a single particle with large momentum, large 
transverse momentum and large IP are used. In total, the sample consists of ~482 000 
signal events. 

Even though the final state is only partially reconstructed due to the missing neutrino, 
the contamination of background is small and the background to signal ratio B/ S is mea- 
sured to be ~ 0.14 in the signal mass region. The main sources of background are events 
containing a D*^ originating from a 6-hadron decay (referred to as D°-from-i?), events 
with a D*~ not from a 6-hadron decay, decays of mesons to the same particles as the 
signal together with an additional pion (referred to as B^) and combinatorial background. 
The different background sources can be disentangled from the signal by exploiting the 
different distributions of the observables m=mKn, ^^=fnKTTn—fnKn, the reconstructed 
B^ decay time t and the mixing state q. The mixing state is determined by comparing 
the flavour of the reconstructed signal B^ at decay time with the flavour indicated by the 
tagging decision (flavour at production time). For unmixed (mixed) events g=+l(— 1) 
while for untagged events g=0. The decay time is calculated using the measured B^ de- 
cay length, the reconstructed B^ momentum and a correction for the missing neutrino 
determined from simulation. It is parametrized as a function of the reconstructed B^ 
invariant mass. 

An extended unbinned maximum likelihood fit is performed by defining a PDF for the 
observables (m. Am, t, q) as a product of one PDF for the masses and one for the t and 
q observables. For the and D*~ mass peaks two double Gaussian distributions with 
common mean are used, while a parametric function motivated by available phase space 
is used to describe the Am distributions of the D^-from-B, and combinatorial background 



7 



components. The decay time distribution of the signal consists of mixed, unmixed and 
untagged events, and is given by 



^s., N / et,^ a{t) {e-'/^BO [I + q^i - 2uj) cosiAnidt)] ^ R{t - t')} if g = ±1 
[t, gj oc I _ ^^^^^^^^^ ^ _ ^ ^ 



(7) 



where Am^ and r^o are the B^-B^ mixing frequency and B^ hfetime. The decay time 
acceptance function is denoted by a{t) and R{t—t') is the resolution model, both extracted 
from simulation. A double Gaussian distribution with common mean is used for the decay 
time resolution model. In Eq. [7] the tagging parameters are assumed to be the same for 
B and B-mesons. 

The decay time distributions for the B~^ and D'^-irom-B background components are 
taken as exponentials convolved by the resolution model and multiplied by the same 
acceptance function as used for the signal. For the prompt D* and combinatorial back- 
ground. Landau distributions with independent parameters are used. The dependence on 
the mixing observable q is the same as for the signal. The tagging parameters ^tag and 
u of the signal and of each background component are varied independently in the fit, 
except for the B~^ background where they are assumed to be equal to the parameters in 
the signal decay. Figure |2] shows the distributions of the mass and decay time observables 
used in the maximum likelihood fit. The raw asymmetry is defined as 

jYunmix/^\ _ jyinix/^N 

^unmixj^^^ _j_ jymixj^^^ ^ ' 

where A^™™ ^jyunmix-j jg ^.j^g number of tagged events which have (not) oscillated at decay 
time t. From Eq. [7] it follows that the asymmetry for signal is given by 

A{t) = (1 - 2uj) cos(Amrft). (9) 

Figure [3] shows the raw asymmetry for the subset of events in the signal mass region 
that are tagged with the OS tagger combination. At small decay times the asymmetry 
decreases due to the contribution of background events, ^ ~ 0. The value of Am^ was 
fixed to Arrid = 0.507 ps""*^ \9\. Letting the Am^ parameter vary in the fit gives consistent 
results. 

4.3 Analysis of the — )- J/t/jK*^ channel 

The 5° — )■ J/ip K*^ channel is used to extract the mistag rate through a fit of the flavour 
oscillation of the B^ mesons as a function of the decay time. The flavour of the B^ meson 
at production time is determined from the tagging algorithms, while the flavour at the 
decay time is determined from the K*^ flavour, which is in turn defined by the kaon 
charge. 

The 5° — )■ J/ipK*^ candidates are selected from J/tp — )■ fi'^fi~ and K*'^ A'+vr" 
decays. The J/ip mesons are selected by the same selection as used for the 5+ — )• J/ijj AT"*" 



channel, described in Sect. 4.1 The K*^ candidates are reconstructed from two good 
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B decay time (ps) 

Figure 2: Distributions of (a) K~^tt~ invariant mass, (b) mass difference m{KTnr)—m(K7r) and 
(c) decay time of the D*^ fi'^v^ events. Black points witli errors are data, the blue curve 

is the fit result. The other lines represent signal (red dot-dashed), D^-hovn-B decay background 
(gray dashed), background (green short dashed), D* prompt background (magenta solid). 
The combinatorial background is the magenta filled area. 

quality charged tracks identified as and tt". The reconstructed K*^ meson is required 
to have a transverse momentum higher than 1 GeV/ c, a good quahty vertex and an in- 
variant mass within ± 70 MeV/ of the nominal K*^ mass. Combinations of J /if) and K*^ 
candidates are accepted as candidates if they form a common vertex with good quality 
and an invariant mass in the range 5100 — 5450 MeV/c^. The transverse momentum is 
required to be higher than 2 GeV/ c. The decay time and the invariant mass of the B^ are 
extracted from a vertex fit with an identical procedure as for the B^ ^ J/ipK^ channel, 
by applying a constraint to the associated primary vertex, and a constraint to the J/ip 
mass. In case of multiple B candidates per event, only the candidate with the smallest 
of the vertex is kept. 

Only events that were triggered by the "lifetime unbiased" selection are kept. The 
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Figure 3: Raw mixing asymmetry of ^ D*^ ^^v^ events in the signal mass region when 
using the combination of all OS taggers. Black points are data and the red solid line is the result 
of the fit. The lower plot shows the pulls of the residuals with respect to the fit. 

candidates are required to have a decay time higher than 0.3 ps to remove the large 
combinatorial background due to prompt J/ifj production. The sample contains ~ 33 000 
signal events. 

The decay time distribution of signal events is parametrized as in Eq. [7[ without the 
acceptance correction. The background contribution, with a background to signal ratio 
B/S ~ 0.29, is due to misreconstructed 6-hadron decays, where a dependence on the 
decay time is expected (labelled "long-lived" background). We distinguish two long-lived 
components. The first corresponds to events where one or more of the four tracks originate 
from a long-lived particle decay, but where the flavour of the reconstructed K*^ is not 
correlated with a true 6-hadron. Its decay time distribution is therefore modelled by a 
decreasing exponential. In the second long-lived background component, one of the tracks 
used to build the K*^ originated from the primary vertex, hence the correlation between 
the K*^ and the B flavour is partially lost. Its decay time distribution is more "signal- 
like", i.e. it is a decreasing exponential with an oscillation term, but with different mistag 
fraction and lifetime, left as free parameters in the fit. 

The signal and background decay time distributions are convolved with the same 
resolution function, extracted from data. The mass distributions, shown in Fig. [4], are 
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described by a double Gaussian distribution peaking at the mass for the signal com- 
ponent, and by an exponential with the same exponent for both long-lived backgrounds. 

The OS mistag fraction is extracted from a fit to all tagged data, with the values 
for the lifetime and Am^ fixed to the world average [oj . Figure [s] shows the time- 
dependent mixing asymmetry in the signal mass region, obtained using the information 
of the OS tag decision. Letting the Amd parameter vary in the fit gives consistent results. 




Figure 4: Mass distribution of OS tagged — )- J/tpK*^ events. Black points are data, the 
solid blue line, red dotted line and green area are the overall fit, the signal and the background 
components, respectively. 



5 Calibration of the mistag probability on data 

For each individual tagger and for the combination of taggers, the calculated mistag 
probability {r}) is obtained on an event-by-event basis from the neural network output. 
The values are calibrated in a fit using the measured mistag fraction (u) from the self- 
tagged control channel B^ — )■ J/ipK^ . A linear dependence between the measured and 
the calculated mistag probability for signal events is used, as suggested by the data 
distribution, 

^{v) = Po + Pi{v - {v)) , (10) 

where po and pi are parameters of the fit and {rj) is the mean calculated mistag probability. 
This parametrization is chosen to minimize the correlation between the two parameters. 
Deviations from po = (^7) and pi = 1 would indicate that the calculated mistag probability 
should be corrected. 

In order to extract the and pi calibration parameters, an unbinned maximum like- 
lihood fit to the mass, tagging decision and mistag probability r] observable is performed. 
The fit parametrization takes into account the probability density function of r/, V{ri), 
that is extracted from data for signal and background separately, using events in different 
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Figure 5: Raw mixing asymmetry of the — )• J/iIjK*^ events in the signal mass region, for 
all OS tagged events. Black points are data and the red solid line is the result of the fit. The 
lower plot shows the pulls of the residuals with respect to the fit. 

mass regions. For example, the PDF for signal events from Eq. |6]then becomes 

{f^tag (1 — ^iji)) 'P^iv) "right tag decision" 
^ta.g ^{v)'P^{v) r= "wrong tag decision" (11) 

1 ~ ^tag r="no tag decision". 

The measured mistag fraction of the background is assumed to be independent from the 
calculated mistag probability, as confirmed by the distribution of background events. 

The calibration is performed on part of the data sample in a two-step procedure. 
Each tagger is first calibrated individually. The results show that, for each single tagger, 
only a minor adjustment of po with respect to the starting calibration of the neural 
network, performed on simulated events, is required. In particular, the largest correction 
is po— {rj) = 0.033±0.005 in the case of the vertex charge tagger, while the deviations from 
unity of the pi parameter are about 10%, similar to the size of the corresponding statistical 
errors. In a second step the calibrated mistag probabilities are combined and finally the 
combined mistag probability is calibrated. This last step is necessary to correct for the 
small underestimation (po — (^7) = 0.022±0.003) of the combined mistag probability due 
to the correlation among taggers neglected in the combination procedure. The calibrated 
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mistag is referred to as rjc in the following. 

Figure |6] shows the distribution of the mistag probability for each tagger and for their 
combination, as obtained for — t- J/ip events selected in a ±24 MeV/ mass window 
around the mass. 




0.1 0.2 0.3 0.4 0.5 0.6 



'c 

Figure 6: Distribution of the calibrated mistag probability for the single OS taggers and their 
combination for S+ — )• J/ipK~^ events selected in a ±24 MeV/c^ mass window around the 
mass. 



6 Tagging performance 

The tagging performances of the single taggers and of the OS combination measured 
after the calibration of the mistag probability are shown in Tables |2| |3] and |4] for the 
J/ipK'^, 5°-)- J/ipK*^ and D*~^^v^ channels, respectively. 

The performance of the OS combination is evaluated in different ways. First the av- 
erage performance of the OS combination is calculated, giving the same weight to each 
event. In this case, the best tagging power is obtained by rejecting the events with a poor 
predicted mistag probability rjc (larger than 0.42), despite a lower Stag- Additionally, to 
better exploit the tagging information, the tagging performance is determined on inde- 
pendent samples obtained by binning the data in bins of rj^. The fits described in the 
previous sections are repeated for each sub-sample, after which the tagging performances 
are determined. As the samples are independent, the tagging efficiencies and the tagging 
powers are summed and subsequently the effective mistag is extracted. The total tagging 
power increases by about 30% with respect to the average value, as shown in the last line 
of Tables Ha 

The measured tagging performance is similar among the three channels. The differ- 
ences between the B+ ^ J/i)K+ and J/ijjK*° resuhs are large in absolute values, 
but still compatible given the large statistical uncertainties of the B^ — )■ J/ip K*^ results. 
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Differences between the tagging efficiency in the — )■ D*^fi^h'^ and the B — )■ J/ipX 
channels were shown in previous MC studies to be related to the different B momentum 
spectra and to different contributions to the trigger decision [sj. 

Table 2: Tagging performance in the S^— )• J/ij)K~^ channel. Uncertainties are statistical only. 



Taggers 




u [%] 




£tag(l - 2coy [%] 




4.8±0.1 


29.9±0.7 


0.77±0.07 


e 


2.2±0.1 


33.2±1.1 


0.25±0.04 


K 


11.6±0.1 


38.3±0.5 


0.63±0.06 


Qvtx 


15.1±0.1 


40.0±0.4 


0.60±0.06 


OS average {ijc <0.42) 


17.8±0.1 


34.6±0.4 


1.69±0.10 


OS sum of rjc bins 


27.3±0.2 


36.2±0.5 


2.07±0.11 



Table 3: Tagging performance in the t- J/i/j K*^ channel. Uncertainties are statistical only. 



Taggers 




u 


[%] 




£tag(l - 2coy [%] 




4.8±0.1 


34.3±1.9 


0.48±0.12 


e 


2.2±0.1 


32.4±2.8 


0.27±0.10 


K 


11.4±0.2 


39.6±1.2 


0.49±0.13 


Qvtx 


14.9±0.2 


41.7±1.1 


0.41±0.11 


OS average (r^^ <0.42) 


17.9±0.2 


36.8 


±1.0 


1.24±0.20 


OS sum of ?7c bins 


27.1±0.3 


38.0±0.9 


1.57±0.22 



Table 4: Tagging performance in the — )• D* fi^u^ channel. Uncertainties are statistical 
only. 



Taggers 


^tag[%] 


u 


[%] 




£tag(l - 2C0Y [%] 


/i 


6.08±0.04 


33.3±0.4 


0.68±0.04 


e 


2.49±0.02 


34.3±0.7 


0.25±0.02 


K 


13.36±0.05 


38.3±0.3 


0.74±0.04 


Qvtx 


16.53±0.06 


41.5±0.3 


0.48±0.03 


OS average (r^c <0.42) 


20.56±0.06 


36.1±0.3 


1.58±0.06 


OS sum of rjc bins 


30.48±0.08 


37.0±0.3 


2.06±0.06 
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7 Systematic uncertainties 



The systematic uncertainties on the cahbration parameters po and pi are studied by 
repeating the cahbration procedure on — )■ J/ip events for different conditions. The 
difference is evaluated between the value of the fitted parameter and the reference value, 
and is reported in the first row of Table |5} Several checks are performed of which the 
most relevant are reported in Table |6] and are described below: 



Table 5: Fit values and correlations of the OS combined mistag calibration parameters measured 
in the — ;> J/ipK^ , — )• J/ipK*^ and B^ — )• Z)*~/i+f^ channels. The uncertainties are 
statistical only. 



Channel 


Po 


Pi 


iVc) 


Po-Pi{Vc) 


p{Po,Pi) 




0.384 ± 0.003 


1.037 ±0.038 


0.379 


-0.009 ±0.014 


0.14 




0.399 ± 0.008 


1.016 ±0.102 


0.378 


0.015 ±0.039 


0.05 




0.395 ±0.002 


1.022 ±0.026 


0.375 


0.008 ±0.010 


0.14 



Table 6: Systematic uncertainties on the calibration parameters po and pi obtained with )• 
J/ipK^ events. 



Systematic effect 


Spo 


6pi 


S{Po-Pi{Vc)) 


Run period 


±0.005 


±0.003 


±0.004 


B-flavour 


±0.008 


±0.067 


±0.020 


Fit model assumptions V{r]) 


< ±0.001 


±0.005 


±0.002 


Total 


±0.009 


±0.07 


±0.02 



• The data sample is split according to the run periods and to the magnet polarity, 
in order to check whether possible asymmetries of the detector efficiency, or of the 
alignment accuracy, or variations in the data-taking conditions introduce a difference 
in the tagging calibration. 

• The data sample is split according to the signal flavour, as determined by the recon- 
structed final state. In fact, the calibration of the mistag probability for different B 
flavours might be different due to the different particle/antiparticle interaction with 
matter or possible detector asymmetries. In this case a systematic uncertainty has 
to be considered, unless the difference is explicitly taken into account when fitting 
for CP asymmetries. 

• The distribution of the mistag probability in the fit model, V{rj), is varied either by 
assuming the signal and background distributions to be equal or by swapping them. 
In this way possible uncertainties related to the fit model are considered. 
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In addition, the stability of the cahbration parameters is verified for different bins of 
transverse momentum of the signal B. 

The largest systematic uncertainty in Table |6] originates from the dependence on the 
signal flavour. As a cross check this dependence is also measured with — ^ D*^/i+z/^ 
events, repeating the calibration after splitting the sample according to the signal decay 
flavour. The differences in this case are 6pQ = ±0.009 and 6pi = ±0.009, where the latter 
is smaller than in the B~^ — )■ J/ip channel. Both for the run period dependence and for 
the signal flavour the variations of 5pQ and 5pi are not statistically significant. However, as 
a conservative estimate of the total systematic uncertainty on the calibration parameters, 
all the contributions in Table |6] are summed in quadrature. The tagging efficiencies do 
not depend on the initial fiavour of the signal B. In the case of the B^ — )■ J/ipK^ channel 
the values are (27.4 ± 0.2)% for the B+ and (27.1 ± 0.2)% for the 5". 



8 Comparison of decay channels 

The dependence of the calibration of the OS mistag probability on the decay channel is 
studied. The values of po, pi and (rjc) measured on the whole data sample for all the three 
channels separately, are shown in Table |5j The parameters pi are compatible with 1, 
within the statistical uncertainty. The differences po — pi(?7c), shown in the fifth column, 
are compatible with zero, as expected. In the last column the correlation coefficients are 
shown. 

To extract the calibration parameters in the B^ — )■ J/ipK*^ channel an unbinned 
maximum likelihood fit to mass, time and rjc is performed. In analogy to the B^ — )■ J/ijj 
channel, the fit uses the probability density functions of rjc-, extracted from data for signal 



and background separately by using the sPlot 10 technique. The results confirm the 
calibration performed in the B^ — t- J/ip channel, albeit with large uncertainties. The 
results for the B^ — )■ D*~fi^u^ channel are obtained from a fit to independent samples 
corresponding to different ranges of the calculated mistag probability as shown in Fig. [7j 
The trigger and offline selections, as well as signal spectra, differ for this decay channel 
with respect to the channels containing a J/V' meson. Therefore the agreement in the 
resulting parameters is a validation of the calibration and its applicability to B decays 
with different topologies. In Fig. |8]the dependency of the measured OS mistag fraction as 
a function of the mistag probability is shown for the B~^ — )■ J/ipK~^ and B^ — )• D*~ ^'^v^ 



signal events. The superimposed linear fit corresponds to the parametrization of Eq. [TO 
and the parameters of Table |5| 

The output of the calibrated fiavour tagging algorithms will be used in a large variety of 
time- dependent asymmetry measurements, involving different B decay channels. Figure [9] 
shows the calculated mistag distributions in the B^ — )■ J/ipK'^, B^ — )■ J/ipK*^ and 
Bg — )• J/ip channels. These events are tagged, triggered by the "lifetime unbiased" lines 
and have an imposed cut of t > 0.3 ps. The event selection for the decay -8°—)- J/'ipcp is 
described elsewhere [3]. The distributions of the calculated OS mistag fractions are similar 
among the channels and the average does not depend on the pi of the B. It has been also 
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Table 7: Tagging efficiency, mistag probability and tagging power calculated from event- 
by-event probabilities for 5+^ J/i)K+, 5°^ J/i)K*^ , B^^ D*-yi+v^ and 5°^ J/V^0 
signal events. The quoted uncertainties are obtained propagating the statistical (ffist) 
and systematic (second) uncertainties on the calibration parameters determined from the 



Channel 




[%] 






[%] 




e,.,V^ [%] 


J/iIjK+ 


27.3 ±0.1 
27.3 ±0.3 
30.1 ±0.1 
24.9 ±0.5 


36.1 ±0.3 ±0.8 

36.2 ±0.3 ±0.8 
35.5 ±0.3 ±0.8 
36.1 ±0.3 ±0.8 


2.10 ±0.08 ±0.24 
2.09 ±0.09 ±0.24 
2.53 ±0.10 ±0.27 
1.91 ±0.08 ±0.22 



checked that the mistag probability does not depend on the signal B pseudorapidity. 



9 Event-by-event results 

In order to fully exploit the tagging information in the CP asymmetry measurements, 
the event-by-event mistag probability is used to weight the events accordingly. The ef- 
fective efficiency is calculated by summing the mistag probabilities on all signal events 
Y,i (1 - 2uo{7]if)/N. We underline that the use of the per-event mistag probability allows 
the effective efficiency to be calculated on any set of selected events, also for non flavour- 
specific channels. Table [7] reports the event-by-event tagging power obtained using the 
calibration parameters determined with the B^ — )■ J/il) events as reported in Table [s] 
The uncertainties are obtained by propagating the statistical and systematic uncertainties 
of the calibration parameters. In addition to the values for the three control channels the 
result obtained for i?^ — )■ J/ip (p events is shown. For all channels the signal is extracted 
using the sPlot technique. The results for the tagging power are compatible among the 
channels containing a J/ip meson. The higher value for 5°—)- D*~ii^u^ is related to the 
higher tagging efficiency. 



10 Summary 

Flavour tagging algorithms were developed for the measurement of time-dependent asym- 
metries at the LHCb experiment. The opposite-side algorithms rely on the pair production 
of h and h quarks and infer the flavour of the signal B meson from the identification of the 
flavour of the other h hadron. They use the charge of the lepton (/i, e) from semileptonic 
B decays, the charge of the kaon from the 6 — )■ c — )■ s decay chain or the charge of the 
inclusive secondary vertex reconstructed from 6-hadron decay products. The decision of 
each tagger and the probability of the decision to be incorrect are combined into a sin- 
gle opposite side decision and mistag probability. The use of the event-by-event mistag 
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probability fully exploits the tagging information and estimates the tagging power also in 
non flavour-specific decay channels. 

The performance of the flavour tagging algorithms were measured on data using three 
flavour-specific decay modes J/iJjK~^, — )■ J/ipK*'^ and — > D*~^~^h'^ . The 

B~^ — )■ J/ip channel was used to optimize the tagging power and to calibrate the mistag 
probability. The calibration parameters measured in the three channels are compatible 
within two standard deviations. 

By using the cahbration parameters determined from S+ — > J/z/; K'^ events the OS tag- 
ging power was determined to be etag(l -2a;)^ = (2.10±0.08±0.24)% in the B+-> J/ipK^ 
channel, (2.09±0.09±0.24)% in the 5°^ J/ijK*° channel and (2.53±0.10±0.27)% in the 
B^ D*~ fi^v^ channel, where the first uncertainty is statistical and the second is sys- 
tematic. The evaluation of the systematic uncertainty is currently hmited by the size of 
the available data sample. 
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B decay time (ps) 



B decay time (ps) 



Figure 7: Raw mixing asymmetry as a function of B decay time in ^ D*~ ii^v^j, events, in 
the signal mass region, using the OS tagger. Events are split into seven samples of decreasing 
mistag probability r/c- 
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Figure 8: Measured mistag fraction (cj) versus calculated mistag probability (r^c) calibrated on 
^ J/tpK'^ signal events for the OS tagger, in background subtracted events. Left and right 
plots correspond to — >• J/ijjK~^ and ^ D*~^'^v^ signal events. Points with errors are 
data, the red lines represent the result of the mistag calibration, corresponding to the parameters 
of Tables 
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